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Abstract 

I Systems built out of N-body interactions, beyond 2-body interactions, are formulated on the plane, 

' and investigated classically and quantum mechanically (in phase space). Their Wigner Functions — the 



density matrices in phase-space quantization — are given and analyzed. 



- T— I I 

^ ■ In this brief note, we consider systems of N particles on a plane, interacting through N-body potentials. 



over and above more familiar 2-body potentials. Below, we provide the compact generating function. 



G a,a ,b,b ;z,z ,p,p = ^^exp- mz-a • 2n(z - a) , (1 



for their Wigner Functions (WF), the density matrices in phase-space quantization [0, ^, ^, (For reviews, 



see 



Consider the hamiltonian of three particles on a plane, interacting in proportion to the (signed) area 
of the triangle they define: 

H, = ii±^i±Ei^„,A, (2) 

2m 

where the signed area (doubled) is 

^=(ri-r2)A(r2-r3). (3) 
^ On leave from the Theoretical Physics Department, Uppsala University, Sweden. 
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IS 



Its (saddle) quadratic form, with the three 2-vectors arrayed in succession, R = (xi, yi, X2, 2/2, a^s, ys); 



(4) 
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This may be simpUfied by ehmination of the 2x2 structure. On the plane, rotational invariance can 
be exploited by introducing complex variables, z = x + iy, which cuts down the size of such matrices by 
half on each side, since 



/ * *^ / 1 
ri Ar2 = a:i2/2 -a;2yi = Y (^I'^a) ( q 



Z2 



(5) 



Clearly multiplication of the zs by an arbitrary phase yields a rotated configuration with the same eigenvalue 
(a double degeneracy of the original problem). 



Thus the above quadratic form hermitean matrix (up to 1/2) reduces to the imaginary antisymmetric 
matrix 



ft 



Consequently, the above hamiltonian amounts to 
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1 * "T- * r-. 

^ 2m ^ ^ 2 

for complex 3-vectors z and their canonical conjugates p* (with components p = Px + ipy) 



(6) 



(7) 



The matrix fl consists of Sylvester's celebrated "nonions" Q, today's standard clock and shift matrix 
basis. Specifically, Q = —i{h — h"^), where h is the cyclic permutation shift matrix, with /i^ = 1, hence 
/i* = = h? . For UJ = — 1/2 + i\/3/2, a cube root of unity, is diagonalized by the Finite Fourier 
Transform unitary matrix ||^], 



U = 

u^m = 

i.e., the (real) difference of two clock matrices. 




(8) 
(9) 



The normal mode frequencies-squared of the hamiltonian then are 0, \/3, and — -v/3. The negative 
eigenvalues (imaginary frequencies) reflect the cyclic symmetry, 1,2,3, anti-mirrored into the antistandard 
order, 1,3,2. For every eigenvalue, the opposite eigenvalue follows for the two-particle permuted (here, 
antistandard) configuration. 



The corresponding eigenvectors describe the translation zero mode, (1,1,1), with all three particles 
moving in coincidence, so without a force; the pulsating cube roots of unity configuration in standard. 



Curtright, Polychronakos & Zachos 



hep-th/Ollim 



Area Potentials 



3 




1,2,3, order, {uj,lj'^, 1), , or any rotation of it on the plane (with which it can be 

combined to rotating configurations); and, finaUy, its mirror image (c.c), (a;^,u;, 1), in antistandard, 1,3,2, 
order. These last (imaginary frequency) unstable modes lead to the indefinite growing of the area of the 
antistandard-order triangle, corresponding to a quantum mechanical spectrum unbounded below. 



The unstable normal modes of increasing negative area and the concomitant lower unboundedness of 
the spectrum may be counteracted by considering additional harmonic two-body interactions between each 
pair of particles, i.e., a positive semi-definite potential component, 

V2 = ^ ((ri - r^f + (r2 - r^f + (rg - r,f) , (10) 
whose corresponding quadratic form matrix, 

projects out the translation mode and has eigenvalue \/3 for the other modes. Thus, the potential A + V2, 
corresponds to the modified form W = + W, whose eigenvalues are thus shifted from the above ones 
to non-negative ones, 0, 2^/3, and 0, for the same orthogonal eigenvectors. Consequently, considering the 
plane-rotational invariance doubling of the modes, the system ^3 resolves to two harmonic oscillators and 
four free modes, whose quantization yields a 2-d linear spectrum on the free particles' continuum. 



The system considered so far naturally generalizes to a standard (anticlockwise ordered) A^-gon, with 
(twice) the area, consisting of a fan of triangles, 

A = ri Ar2 r2 A rs ... TAT Ari. (12) 

It corresponds to the N x N matrix ^, ^ 

n = -i{h - h^-^).... = 2 sin (ln(-i/i)) , (13) 

where h is the N x N cyclic shift matrix, with = 1 and = h^^ = h^^^ . Extension from to -|- 1 
amounts to the addition of just one triangle to the fan. 



Thus, the general eigenvalues for $7 are 

Afc = -i(u;'=-0, (14) 

with oj the primary A^-th root of unity and A; = 0,1,2,...,A'^ — 1. Note the evident mirror-image pairing of 
the eigenvalues = —\N-k- The corresponding N complex eigenvectors are 

|efc) = (l,a;^a;2^..., 0.^(^-1)), (15) 

which array into the corresponding unitary diagonalizing Finite Fourier Transform matrix U . 
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i<j \ i i / 



The negative eigenvalues may be shifted, as above, by the addition of a 2-body interaction in the 
potential, 

^2 = 1 - = ^\Ny:r^- (> Jr,)^ ) , (16) 

i<j 

where c is a normahzation constant. This interaction corresponds to the matrix 

W = -{N l-|eo)(eo|). (17) 
c 

I Co) = (1,1,1,1,...,1) is the translation zero mode projected out, while all other eigenvalues but its own 
are shifted by N/c. Thus, choosing the normalization (c < N/2 will suffice), all negative eigenvalues of 213 
may be neutralized to yield a stable system Sjn with a spectrum bounded below. 

The classical equations of motion are 

z = -nz. (18) 

Note the vanishing of the torque (and hence conservation of angular momentum) for the system Hi\f 
(manifestly rotational invariant!), even though the interaction is not two-body. This follows from summing 
the torque on each particle, since ft is hermitean, 

r = ^ ri A mi^i = — (z* • fiz - (fiz)* • z) = 0. (19) 



Further note the N independent time- invariant complex vector combinations, 

ex.p{-iVn t) (p + z) . (20) 



In the less compact notation of 2n-dimensional phase space (n = 2A^) there are 2n real independent time- 
invariant combinations. Eliminating t among them yields 2n — 1 independent real time-invariant quantified 
Qi, which characterize maximally superintegrable systems in phase space [^. For these, time evolution is 
completely specified by flow perpendicular to all the phase-space gradients VQj. Thus, the time evolution 
of the phase-space variables, and so the time derivative of any phase-space function, is proportional to the 
corresponding Jacobian determinant for the full 2n-dimensional phase space, 

df d{f,Qi,...,...,Q2n-i) 

dt d{xi,Pa:i,yi,Py^,...,XN,PxM^yN,PyN) 

i.e., the phase-space Nambu Bracket 

There is an alternate way to stabilize the system by shifting the negative eigenvalues of Q without 
the above introduction of 2-body interactions, but, instead, through the introduction of a magnetic field 
B pointing into the plane in question, and consideration of the particles as charged. Such a magnetic 
field breaks parity and differentiates between left- and right-spinning cyclotron orbits. The magnetic field 
modifies the complex equations of motion ( |T8[ ) , 

i = iB±- ftz, (22) 



^ For example, in the normal mode representation of f)3 with the full six-dimensional 2IJ diagonalized to 
diag(0, 0, 2\/3, 2\/3, 0, 0), and the coordinates starting at at f = 0, a set of 11 independent invariants is 

p^^ + »2%/3a;2 ^ {p +i2V3x2)exp{~2V3xi/mp^^)- Px^, Py^, Px^, Pys', Evilly P^afi^ Pm^ 
Py^+i2VSy2 PxiVi Px^xs p^j^ya 
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and thus shifts the eigenfrequencies of the modes from -/Xk to Wk, solutions of 

wl - Bwk - Afc = 0. (23) 

The modified mode frequencies 



Wk = ^±]J— + Xk (24) 

are real for B^ > — 4min(Afc), which equals 8 for even A'^, and is smaller for odd N. Thus, the system 
is stable for a field at least as strong as \B\ > 2\/2. (In the small mass limit, where the strength of the 
area potential and the magnetic field are kept unchanged, the left-hand side of (22) drops out and hence 
Wk = ^k/B^ describing pure Landau levels.) 



The zero mode = coincides with yAo = and represents Landau level degeneracy. The mode 
frequencies Wk can be positive or negative, corresponding to the chirality of the cyclotron orbits. At the 
critical field, \Bc\ = 2\/2, (taking —iuik = ex'p{27rik/N — i7r/2)), the above square root acts on a perfect 
square, so that 

Wk/2 = sgn(5,) cos(^/4) ± cos - , (25) 

where sgn(i?c) = ±1. So, essentially (for large where the vanishing and nonvanishing w corresponding 
to Aq = are neglected), 25% of the modes are of one chirality and 75% are of the opposite. For the above 
= 3 case, the critical field is Be = 2(3)^/^; beyond the zero mode, there is one mode with w = Be, two 
with w = (3)^/'^(l ± and the two ones stabilized by the field to it; = Be/ 2. 



In the large A'^ limit, the above system reduces to a closed noncovariant string undulating on the plane. 
The particle index becomes a continuous periodic variable, a = 27rn/N, and the complex A^- vector of the 
particles' position goes into a scalar closed string field of that variable, z(t) i-^ (p{a, t). The potential is the 
(signed) area enclosed by the string. The resulting noncovariant lagrangian density is 

tn — • — 

C = -<i)(t> + i(t)d„<i). (26) 

This amounts to Schrodinger's Lagrangian, with the roles of space and time reversed. The dispersion 
relation of Schrodinger's action, E = p^/2, here reverses to 

w'^ + 2k = 0, (27) 

with k the integer Fourier modes of the compact variable a. This specifies an infinite tower of stable 
and unstable modes. As in the finite case, introducing an additional harmonic two-body potential for the 
particles amounts to a term, 

2 

(28) 

essentially an external harmonic potential minus a center-of-mass term. This would correspond to an 
external constant potential in the Schrodinger equation. The modes are shifted by a constant, but still an 
infinity of them remain unstable. 



<T<cr' 



d(jd(j'\cp{a) - (/.(ct')P = j d(j\^{a)\'^ - J da(f){a 



We finally quantize the original system in phase-space [^. Groenewold's associative ★-product |T^, in 
our condensed notation of complex vectors. 



:,ih{dz*-dp+dz-dp* — dp*-dz.— dp-dz*) 



(29) 
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is the cornerstone of phase-space quantization. (It allows c-number functions in phase space to multiply 
with each other and with WFs associatively and noncommutatively, in perfect parallel to operator ma- 
nipulations in the standard, Hilbert space, formulation of quantum mechanics ||5|.) The Poisson Brackets 
which are exponentiated in the ^-product are 

{z*,p,} = 26ij , (30) 

and their complex conjugate, for i,j = 1, N. 

The linear contact transformation U discussed above, which leads to the normal mode variables U^z 
and f/tp, is a linear canonical transformation, and therefore preserves Poisson Brackets and ^-products 



|11|. (To be contrasted to nonlinear canonical transformations which actually transform ^-products in 
a covariant fashion |]l2|). That is to say, for these transformations the * above is a scalar, i.e. equals 
to the same expression for the normal mode variables U^z. All formal manipulations below, then, may 
be equivalently conducted in the space of the original or the normal mode variables, and no distinction 
between either variables in * is necessary. 

Thus, e.g., in terms of the normal modes, the WF for the = 3 case ^3 is a product of the six WFs 
of each normal mode: four unnormalizable free ones and two oscillators of frequency-squared 2^/3 = w'^, 

f = 6{pi-hk,)6{p2-hk2)^-^^ exp f-i(-^ + mwrl)) L„ f ^(^ + mwxl)) f ^(^ + mwyl)] , 



vr^ V ''TT-w J \ii mw J \h mw 

(31) 

where the L^s are Laguerre polynomials with n,s = 0,1,2,...., and the two ks comprise four arbitrary 
real constants. These WFs obey the phase-space stargenvalue equation for the complete spectrum of the 
quantum system |12], 

f = ( h^^^^ + (^X^\{l+n + s)] f. (32) 



2m V 4 



Nevertheless, the generic case may be approached in a compact way through the WF generating 



functions of |13], which generalize the WFs for coherent states ||TI[ (the multiparticle generalization for 
quadratic systems may also be found in fl^ ). 



G(a, a*,b,b*;z,z*,p,p* 
For 

this function satisfies 



TT 



1 l/|a|2 + |b|2 (p-b)*-(p-b) \ 

m(z — a) • 2li(z — a) ) . (33) 



m 



2m^ ^ 2 



(34) 



^ * G ~ - (2Ar + (a - b) • (^a - Sb) + (a* - b*) • (^a- - )) G, 



G*^ ~ - (2Af + (a + b) • (5a + ^b) + (a* + b*) • (^a. + 5b.)) G. 



(35) 
(36) 



The SB may be diagonalized by unitarily transforming all variables and shifts z, p, a, b, dz* , dp* , through 
U^z, and their c.c.s through Uz*, Udz, etc. All quantities Sj,G,-k, are scalar dot products. Applications of 
such generating functions are detailed in p^ . 
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Several constructions based on the potentials introduced here will be presented in future publications. 
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